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ABSTRACT 

We present rigorous simulations of EUV masks with technological imperfections like side- wall angles and corner 
roundings. We perform an optimization of two different geometrical parameters in order to fit the numerical 
results to results obtained from experimental scatterometry measurements. For the numerical simulations we use 
an adaptive finite element approach on irregular meshes. 1 This gives us the opportunity to model geometrical 
structures accurately. Moreover we comment on the use of domain decomposition techniques for EUV mask 
simulations. 2 Geometric mask parameters have a great influence on the diffraction pattern. We show that 
using accurate simulation tools it is possible to deduce the relevant geometrical parameters of EUV masks from 
scatterometry measurements. 

This work results from a collaboration between AMTC (mask fabrication) , Physikalisch-Technische Bundes- 
anstalt (scatterometry) and ZIB/ JCMwave (numerical simulation). 
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1. INTRODUCTION 

Extreme ultraviolet (EUV) lithography is considered as the main candidate for further miniaturization of com- 
puter technology. Since compared to state-of-the art photomasks, EUV masks are illuminated at oblique in- 
cidence, the quality of pattern profiles becomes important due to shadowing effects. 3, 4 Consequently, there 
is a need for adequate destruction free pattern profile metrology techniques, allowing characterization of mask 
features down to a typical size of 100 nm. 

Here we present an indirect method for the determination of geometrical EUV mask parameters. Ex- 
perimental scatterometry measurements are compared to numerical simulations of EUV masks using the finite 
element method (FEM). 
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Figure 1. (a) Experimental setup for scatterometry experiment with fixed incident angle of 6° and variable angle of 
detection out ; (b) Result of single wavelength scatterometry measurement at A = 13.65 nm. Diffraction orders appear as 
peaks with finite width, the zeroth diffraction peak is centered around 6°. 

2. CHARACTERIZATION OF EUV MASKS BY EUV SCATTEROMETRY 

Single wavelength scatterometry, the analysis of light diffracted from a periodic structure, is a well suited tool 
for analysis of the geometry of EUV masks. Since scatterometry only needs a light source and a simple detector 
with no imaging lens system, its setup is inexpensive and offers no additional technical challenges. Fig. 1(a) 
shows a sketch of the experimental setup. Light of fixed wavelength and fixed incident angle is reflected from 
the mask and the intensity of the reflected light is measured in dependence on the diffraction angle. The use 
of EUV light for mask characterization is advantageous because it fits the small feature sizes on EUV masks. 
Diffraction phenomena are minimized, and of course the appropriate wavelength of the resonant structure of the 
underlying multilayer is chosen. Light is not only reflected at the top interface of the mask but all layers in the 
stack contribute to reflection. Therefore one expects that EUV radiation provides much more information on 
relevant EUV mask features than conventional long wavelength methods. 

All measurements for the present work were performed by the Physikalisch-Technische Bundesanstalt (PTB) 
at the electron storage ring BESSY II. 6 PTB's EUV reflectometer installed in the soft X-ray radiometry beamline 
allows high-accuracy measurements of very large samples with diameters up to 550 mm. 7-9 



3. FEM SIMULATION OF EUV SCATTEROMETRY 

Fig. 1(b) shows the result of a scatterometry measurement of an EUV test mask (see table 1) considered in the 
present work. The position of the diffraction angles provide information about the pitch of the EUV absorber 
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Design parameter 


Tcstmask 


Absorber stack sidewall angle a 
Pitch 

Top critical dimension 


90° 
840nm 
140nm 



Table 1. Design parameters (see also Fig. 4) for EUV test mask produced by AMTC. 



pattern. However the intensities of the diffraction orders do not carry direct information about other topological 
features of the mask. The determination of these features from a scatterometry measurement is a challenging 



inverse problem and a hot topic of actual research. Accurate and fast numerical simulation of the scattering 
experiment thereby plays a vital role. The FEM method is particularly suited for this application. It has several 




Figure 2. SEM pictures of EUV mask patterns and corresponding triangulated geometries for FEM computation. 
advantages 10 : 

• Maxwell's equations describing the scattering problem are solved rigorously without approximations. 

• The flexibility of triangulations allows modeling of virtually arbitrary structures, as illustrated in Fig. 2. 

• Adaptive mesh-refinement strategies lead to very accurate results and small computational times which 
are crucial points for application of a numerical method to the characterization of EUV masks. 

• Choosing appropriate localized ansatz functions for the solution of Maxwell's equations physical properties 
of the electric field like discontinuities or singularities can be modeled very accurately and don't give rise 
to numerical problems, see Fig. 3. 

• It is mathematically proven that the FEM approach converges with a fixed convergence rate towards the 
exact solution of Maxwell-type problems for decreasing mesh width of the triangulation. Therefore it is 
easy to check if numerical results can be trusted. 




Figure 3. FEM solution for the electric field propagating through a phase mask. The electric field has singular behaviour 
at corners of the absorber and discontinuities at material interfaces. 

Throughout this paper we use the FEM solver JCMharmony for numerical solution of Maxwell's equations. 
JCMharmony has been successfully applied to a wide range of electromagnetic field computations including 
waveguide structures, 11 DUV phase masks, 10 and other nano-structured materials. 12 ' 13 It provides higher 



order edge elements, multigrid methods, a-postcriori error control, adaptive mesh refinement, etc. Furthermore a 
special domain decomposition algorithm implemented in JCMharmony is utilized for simulation of EUV masks. 2 
Light propagation in the multilayer stack beneath the absorber pattern can be determined analytically. The 
domain decomposition algorithm combines the analytical solution of the multilayer stack with the FEM solution 
of the absorber pattern, dramatically decreasing computational time and increasing accuracy of simulation 
results. 



4. DETERMINATION OF EUV MASK PARAMETERS WITH FEM SIMULATION 
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Figure 4. Parameters of EUV mask pattern: (a) Absorber stack sidewall angle a and top critical dimension topCD; (b) 
absorber edge Radius R 

The idea of characterising EUV masks with scatterometry measurements and FEM simulations is the com- 
parison of the intensities I n of experimental diffraction orders n with numerically obtained results. 5 First 
the geometry of the EUV mask is modeled using a finite number of parameters. Then scattering from the EUV 
mask is simulated and the intensities /^ lm of the diffraction orders are computed. The deviation £ between ex- 
perimental and numerical intensities is computed and the parameters which minimize this deviation determined. 
This leads to a finite dimensional optimization problem. In the experimental setup light of each diffraction order 
is always reflected into a finite solid angle leading to peaks with finite width, see Fig. 1(b). The experimental 
intensities 7n Xp that are used are given as the heights of these peaks which are only proportional to the whole 
intensity diffracted into an order (i.e. the integral over a peak). Numerically we determine the whole intensity 
of each diffraction order and therefore we have to scale the simulated intensities /® lm uniformly with a factor 7 
before determining the deviation £: 

where the global scaling factor 7 is determined by minimizing £ 2 with respect to 7: 
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The described procedure was applied to an EUV test mask produced by AMTC. The unknown mask param- 
eters of interest were the sidewall angle of the absorber stack (which was restricted to a < 90°) and the top 
critical dimension (defined in Fig. 4). 

The search for the optimal parameter set of the layout was performed using the Nelder-Mead simplex al- 
gorithm. As starting point the AMTC design parameters were chosen. In order to evaluate the results of our 
procedure scatterometry measurements and FEM simulations were compared at three different EUV wavelengths, 
shown in Fig. 5. 
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Figure 5. Brightfield measurement of multilayer of AMTC test mask: Reflectivity R of open multilayer in dependence on 
incident wavelength for fixed incident angle of 6° . Determination of EUV mask geometry was performed at three spectral 
wavelengths Ai = 13.4 nm, A2 = 13.65 nm, A3 = 13.9 nm. 
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Figure 6. Dependence of deviation £ 2 (see Eq. 1) on geometrical parameters of EUV mask. Fixed parameters: (a) 
a = 90.0°, R = Onm; (b) top CD=146.5nm, R = Onm; (c) a = 90.0°, top CD=146.5nm. 



A comparison of diffraction orders obtained from scatterometry and FEM simulation is shown in Fig. 7 for 
the optimal EUV mask parameters found during optimization. We see excellent agreement at all wavelengths of 
incident EUV radiation. Only the tenth diffraction order at A3 = 13.9 nm differs. Here the simulated intensity 
is much lower than the experimental and therefore the deviation £ 2 much larger than for Ai = 13.4 nm and 
A2 = 13.65 nm, see table 2. The geometrical parameters which belong to the best fitting simulations are shown 
in table 2 in comparison to the desired design values for the mask. We see that the best fitting geometries agree 



geometrical parameter 


design value 


FEM Ai = 13.4 nm 


FEM A 2 = 13.65 nm 


FEM A 3 = 13.9 nm 


a[°] 


90 


87.9 


90.0 


90.0 


top CD [nm] 


140 


145.7 


146.5 


146.5 


deviation £ 2 




0.18 


0.24 


1.32 



Table 2. Deviation £ 2 and geometrical parameters of EUV mask obtained from FEM computation at different wavelengths. 



extremely well for A2 = 13.65 nm and A3 = 13.9 nm. For Ai = 13.4 nm the absorber angle is 2° (2.2%) smaller 
and the top CD is 0.8 nm (0.5%) smaller than for A2 and A 3 . 

Fig. 6(a), (b) shows how the deviation £ 2 between experimental and simulated diffraction orders depends on 
the top critical dimension and the absorber angle a. It grows strongly with increasing distance from the optimal 
geometrical parameters. This shows that the presented method is very robust. 

As a further geometrical parameter the absorber edge radius R was considered, see Fig. 4(b). The best 
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Figure 7. Comparison between experimental scatterometry measurement and FEM computation of diffraction orders for 
different wavelengths of incident EUV light. 
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Figure 8. Convergence of FEM method: relative error of intensity of first two diffraction orders in dependence on number 
of unknowns of FEM computation. 



fitting value for R was determined at the incident wavelength A2 = 13.65 nm and the optimal values for top 
CD = 146.5 nm and a = 90.0°. Fig. 6(c) shows a minimal deviation between experiment and simulation for 
R = 6.2 nm. We see that the edge radius does not have such a great effect on the diffraction orders since the 
deviation £ 2 hardly changes compared to the effects of the top CD (a) and absorber angle (b). This confirms 
further that for the determination of more sophisticated geometrical parameters very accurate simulations are 
crucial. As already mentioned the convergence of the FEM method is mathematically proven and it is therefore a 
very good choice for the presented method. Fig. 8 shows the convergence of the zeroth and first diffraction order. 
We see the relative error in dependence on the number of unknowns of the FEM computation (i.e. a coarser 
triangulation) . Furthermore the computational time on a standard PC (3.4 GHz Intel Pentium 4, 1GB RAM) is 
shown at each refinement step of the grid. After 72 s we already have a relative error of 10~ 3 much smaller than 
the experimental uncertainty of about 0.01. A short computation time also becomes crucial for the determination 
of mask parameters when choosing a larger number of independent geometrical parameters and performing the 
search for the optimal values in a higher dimensional space. We expect that scatterometry measurements at 
several different wavelengths will become very important for the presented method when characterizing EUV 
masks in greater detail. In order to further validate the geometrical parameters of the EUV mask obtained via 
scatterometry and FEM simulation a comparison to direct measurements like atomic force microscopy is planned. 
These measurements will be carried out at AMTC. 



5. CONCLUSIONS 

We demonstrated that single wavelength scatterometry in combination with FEM simulations is a promising 
candidate for an accurate and robust destruction free characterization of EUV masks. Thereby experimental 
diffraction orders are compared to FEM simulations of EUV masks. For FEM simulations the EUV mask is first 
described with a finite number of geometrical parameters like sidewall angles, line widths, corner roundings, etc. 
and then the best fitting values determined by minimizing the deviation of experimental and numerical data. 

Here we considered the top critical dimension, the sidewall angle and the edge radius of the absorber stack 
of an EUV mask as unknown geometrical parameters. The search for the best fitting geometry at three different 
wavelengths gave nearly the same values for the top critical dimension and the absorber sidewall angle proving 
both robustness and accuracy of the method. Furthermore the absorber edge radius had only minor influence 
on the numerical diffraction pattern. 

We showed that very accurate numerical simulations are crucial for detailed geometric characterization of 
EUV masks using scatterometry data. The FEM method is well suited for the simulation of EUV masks since it 
allows computation of nearly arbitrary geometries, is very accurate and very fast. Thereby very fast simulation 
of the EUV mask with a fixed parameter set provides a precondition for the solution of the given inverse problem. 
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